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The complex system of the CMS all-silicon Tracker, with 15 148 silicon strip and 1440 silicon pixel modules, requires sophisticated 
alignment procedures. In order to achieve an optimal track-parameter resolution, the position and orientation of its modules need to 
be determined with a precision of few micrometers. The alignment of pixels modules is crucial for the analyses requiring a precise 
vertex reconstruction. The aligned geometry is based on the analysis of several million reconstructed tracks recorded during the 
commissioning of the CMS experiment, both with cosmic rays and with the first proton-proton collisions. Statistical precision of 
the alignment of the module with respect to the particle trajectories to less than 10 microns has been achieved. The results have 
been validated by several data-driven studies (track fit self-consistency, track residuals in overlapping module regions, and track 
parameter resolution) and compared with predictions obtained from a detailed detector simulation. 

Recent developments include the determination of sensor bow and displacements between sensors of composite modules. 

Thoughts on improving future detectors with respect to alignment are given. 

i Keywords: Pixel detector, CMS, track-based alignment, Millepede-II, Broken Lines 



,0 

on 
© 

© 

s 



1. Introduction : 

The all-silicon inner tracker of the CMS detector at CERN . 
consists of 15 148 silicon strip and 1440 silicon pixel modules . 
in a barrel-and-endacp configuration!]]. Its main purpose is . 
to determine track parameters of charged particles produced . 
in proton-proton and heavy ion collisions. The parameters to , 
be delivered per track are the charge-signed curvature k (i.e. , 
inverse transverse momentum), the impact parameters in the , 
transverse plane and along the beam axis, d xy , d z respectively, , 
and the polar angles 9 and <p. The intrinsic hit resoultions of , 
the detector modules are of the order of a few tens of microns, 
depending on module type and location. In order to determine 
the track parameters with high precision, the positions of the 
sensor modules need to be known better than their resolution. 
Alignment using large amounts of track data (typically several 
millions) is one approach to fulfill this requirement. 



J> 2. Track-based alignment 
• i— i 

The use of tracks to align a tracking detector is possible un- 

p3 der the assumption, that tracks may be described using a limited ; 

20 and sufficient number of parameters in an appropriate way to : 

21 predict their paths. Misalignment leads to a systematic distor- : 

22 tion of the measurements per module which can be determined . 

23 using sufficiently large number of tracks and their hit signals. 
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In CMS, the alignment software consists of two independent 
algorithms, tools for the study of random and systematic mis- 
alignments and an extensive collection of tools to monitor and 
visualize the performance and geometry of the detector. They 
use track data from collision or cosmic ray muons, both as sim- 
ulated and real data. 

Track-based alignment relies on a suitable description of the 
track and its propagation through matter, as defined by the cho- 
sen detector geometry. It can be formulated as a linear least 
squares problem where the following expression needs to be 
minimized: 



tracks hits 



* 2 (P' q) = ZZ r u ( P' y u ri ;(P' ^ 



(i) 



where r, ; is the residual vector containing all residuals from the 
tracks used and their hits, defined as 

r (/ = track-model prediction - measured hit. 

The residuals are a function of p, the vector containing all align- 
ment parameters describing the actual geometry and qj, the 
track parameters of the j track. Vy is the inverse covariance 
matrix containing all information on the measurement precision 
and their correlations. Position and orientation of the detector 
modules contribute 6 or 5 degrees of freedom for silicon pixel 
and strip detectors, respectively. This defines the size of a sub- 
vector of p describing one sensor. 

Using a sufficiently large sample of tracks, equation ([TJ and 
its summands follow a.y 2 distribution for a corresponding num- 
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47 ber of degrees of freedom (ndof), obeying 91 

92 

95 

48 In the case of a detector of the size as in CMS, alignment al- 96 

49 gorithms need to reduce their complexity while preserving their 97 

50 focus on solving the problem for two resons: 1) The total num- 

51 ber of parameters p and q gets large. Aligning 16 000 mod- 99 

52 ules for position and angle leads to about 60000 parameters in 

53 p. The developments described in this paper will increase this 

54 number. 2) The results should be delivered within a reasonable 

55 time-frame. The two approaches in CMS to manage this are as 

56 follows: 

101 

57 Local algorithm. This algorithm reduces the workload by 102 

58 aligning the modules independently ignoring correlations at 103 

59 first glance. Each module is forced to the position predicted™ 
so by the track hits from the other (unaligned) modules. Intrinsi-™ 

61 cally, it uses the track parametrization and propagation from theme 

62 CMS tracking, which takes all necessary effects of the magnetic™ 

63 field and material interactions into account. This approach re-™ 

64 covers correlations between modules by iterating over the same™ 
es event sample several times. At each iteration the tracks are re-"° 
ee fitted using the alignment corrections obtained in the previous"! 
67 iteration. Eventually the procedure yields converged alignment" 2 
es constants. This algorithm is known as the HIP-algorithm$2^\3\."3 

114 

69 Global algorithm. This approach reduces the complexity by" 5 

70 the observation that it is sufficient to solve for the module pa-" 6 

71 rameters p alone. This can be achieved by requiring indepen-"? 

72 dent measurments and the use of block-matrix theorems. This" 6 

73 is implemented in the Millepede-II alignment algorithm[4|. As"9 

74 this presentation focuses on some recent results produced using™ 

75 this algorithm, a more detailled description follows now. ^ 

122 

76 2.1. Millepede-II 123 

124 

77 To accommodate for nonlinearities introduced by the track 

78 parametrization (q) and by the module parameters (p), equation 

79 ([TJ needs to be linearized: 

2 1 / attj atij \ 2 

r(P.q) = Zj Zj^ K " f yfeo.qjo) - ^-Ap - — Aq, 
j i U v 1 1 ' ' 

(3),25 

so where f (J is the hit position predicted by the track model from 
si track reconstruction and my is the measured hit position. As- 127 

82 suming uncorrected measurements allows to replace the in- 

83 verse covariance matrix by -\ with o~i the Gaussian error of 

J o — . 129 

IJ 

84 the measured hit position. 130 
as The track model used in CMS is the Kalman filter descriptioni3i 
se including proper description of material effects [5 1 and the prop-™ 
87 agation in the magnetic field[6|. By design, it is a sequential fiti33 
as and cannot produce the covariance matrix for all track parame-134 

89 ters. In principle, it is possible to gain this information for alli35 

90 tracks a posteriori, but the Broken Lines approach as describedi36 



in Q can be implemented more efficient and is equivalent to 
the Kalman approach. A brief description follows: 

A charged particle traversing material experiences multiple 
scattering, mainly due to Coulomb interaction with the elec- 
trons in the atoms, resulting in a spatial shift and a change of the 
particle direction after leaving the material compared to propa- 
gation in vacuum. The mean of the deflection angle due to this 
effect is (J3) = 0. The distribution of the deflection angles can 
be approximated within certain limits as a Gaussian standard 
deviation cr(B) by the following formula|8 1: 

_ 13.6 MeV ^ j-j^ fi + Q 03g ^ 
vp 

where v - Be (here B as rel. velocity factor) is the velocity of the 
particle, p its momentum and z the charge. x/Xq is the thickness 
of the traversed medium in units of radiation lengths. 

Equation (|4]) takes into account all material traversed by 
the particle for the full trajectory. Care has to be taken dur- 
ing propagation, as simple summing up contributions of sub- 
paths leads to too large estimates of cr(B) by the ln-term in the 
bracket (details in reference). It is standard procedure to treat 
a "thick" scatterer (material with a finite thickness) as two in- 
finitely "thin" scatterers with same mean and sigma spaced by 
1 / V3 of the length of the "thick" scatterer. In a tracking de- 
tector as in CMS, most of the material is concentrated at layers 
coinciding with the detector modules. They consist of matter 
in which the sensor interaction takes place and of non-sensing 
matter like support structures, cabling and cooling pipes. So the 
two "thin" scatterers coincide in the detector planes. 

To determine the momentum of the charged particles' tracks, 
a strong and sufficiently homogeneous magnetic field of 3.8 T is 
present in the tracker. This can be taken into account by adjust- 
ing the expectation value of the scattering angle of a propagated 
particle (B) (=0 without B-field) to the value defined by the ac- 
cumulated Lorentz force while propagating through the field. 

Taking all this into account, the sum over all hits of one track 
in ([3) becomes 

1=1 

«scat-l 

+ A(*,u) r V^,-(*,u) (5) 

i=2 

where /}, is a vector of additional parameters of the track at 
every scatterer to account for the deflection angles, k is the 
charge-signed curvature, u = (ui, . . . u„ sail ) describe the hit po- 
sition in some local frame of the sensor and the projection ma- 
trix P, translates between the track frame and the local frame. 
The sums run over n meas recorded hits and n scat scatterers along 
one track, normally « mea s < "scat as the detector is neither fully 
hermetic nor efficient. 

Solving for the minimum of eq. |5]) leads to a bordered band 
matrix: One /?, depends on the hit and its neighbours only, lead- 
ing to a band matrix structure of band width m. The border 
b in the matrix comes from k, which is connected to every hit 
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along a track. This structure allows for fast solution and de-iss 
termination of the covariance matrix using root-free Choleskyiso 
decomposition with a numeric complexity of 0(n 2 (m + b), com-ioo 
pared to (9(n 3 ) for inversion. This is needed for the refit internalisi 
to MillePede for single tracks. 

All this leads to a track description equivalent to the Kalman 
filter model, as shown in (5|, It has the advantage, that the in- 
verse covariance matrix for one track is a bordered band matrix, 
which can be inverted by root-free Cholesky decomposition, a 
faster approach than inversion. 

3. More detained surface description 

Millepede-II uses an internal track refit as part of its proce- 
dure. Careful studies of its results as a function of track param- 
eters have been carried out. Deviations were found while inves- 
tigating tracks from cosmic ray muons. A strong dependence 
of the ix 2 } on the distance of closest approach do to the beam- 
line (corresponding to track parameter d xy ) was found. This can 
be seen in figure [2j curve for flat module. Several hypotheses 
for the source of this effect have been analyzed. It is an intrin- 
sic property of cosmic rays to have do » in general, which 
translates to a large incident angle on rectangular modules in 
the barrel-shaped part of the detector. For example, whenever 
do is slightly smaller than the radius of a barrel layer, the track 192 
angle becomes especially large, making the hit position highly 193 
sensitive to deviations from an ideal flat rectangle. In the end- 194 
cap parts of the detector, the incident angle of the cosmic ray 195 
muons is always large, independent of do. 196 

To take deviations from a flat surface into account, the sur- 197 
face has been expanded in terms of two-dimensional Legendre 198 
polynomials: 199 



w(u, v) = 2_j c ij L A u ) L i-j( v ) 



(6) 2 



where w{u, v) is the deviation from a plane at the origin of a 
right-handed local cartesian coordinate system (uvw) in w di- 
rection as function of u, v. N is the maximal order of the Leg-204 
endre polynomials. For N — > 00 every possible surface may 
be described, cy are the orthogonal coefficients and L,(x) the2os 
Legendre polynomial of 2-th order. 206 
N = 1 corresponds to the flat module assumption as used207 
prior to this extension (except for translations in u,v and rota-208 
tions around w, which are not covered by |6]) and need the same2oo 
treatment as before). N — 2 introduces three additional parame-210 
ters per module: C20, cn and coi- The first and last can be trans-211 
lated by choosing a proper normalization to sagittae. cn de-212 
scribes a mix-term similar to a twist. Extending the alignment2i3 
to these parameters lead to the curve for flat module in figure2i4 
[2] Figure 1(a) clearly shows how the residuals as a function of2is 



the hit position along strip modules recover to an expected flat2i6 
distribution when taking the bowing of the sensors into account. 217 
The sensors at larger radii of the barrel and endcaps are com-218 
posite modules, i.e. two sensors are mounted in one module2io 
frame and daisy chained to one readout electronic block. Do-220 



shows another deviation from the expected curve. This is at- 
tributed to angles between the sensors. This has been imple- 
mented as well by treating the individual sensors as separate 
entities, which lead to the remaining curves in figure [2] 
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ing the same graph for composite modules (figure 1(b) 1 clearly22i 
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Figure 2: Distribution of the probability of the \ 2 vs - do (MillePede-II track 
refit): For each track the probability of the^ 2 of the MillePede internal track 
refit for the given number of degrees of freedom is calculated. This plot shows 
the average per bin, the error bars are the error of the average. The binning is 
done via the closest distance do of the track to the nominal beamline. Data: 
Cosmic ray muons, recorded in 2009 during commissioning of CMS. 200 000 
tracks used. The results are shown for several cases with different levels of 
description (see text). At large do S 50 cm other effects from the track recon- 
struction start to dominate, which is beyond the normal use-case for the tracker 
and therefore neglected. 



Solving the alignment problem for all these added parame- 
ters lead to the determination of roughly 200000 parameters 
for the full tracker with Millepede-II in one run. This was per- 
formed on a computer equipped with a Intel Nehalem processor 
and 24 GB of RAM within 6 hours of wall-clock time. Crucial 
parts of the algorithm were rewritten for multi-threading using 
OpenMP™ IfTTll to benefit from parallel processing on the 8 
cores the processor offers. The memory consumption for stor- 
ing the matrix of the normal equations was reduced by using 
sparse matrix storage schemes and adaptive selection of storage 
precision of the floating-point numbers at runtime, preserving 
the required overall precision. 

3.1. Estimation of parameter precision 

The Gaussian error of the parameters for the bows were esti- 
mated using the following observation: When solving for a lin- 
ear least squares problem on a computer, the crucial step takes 
place while solving for x in Mx = y, M being the Jacobian 
matrix of the normal equations. M -1 would be the covariance 
matrix of the parameters, usually not feasible to solve for as 
the numerical complexity goes with 0(n 3 ) for matrix inversion, 
compared to other methods for solving for x. For this reason, 
MillePede uses the MINRES algorithm [ 1 1 as a solver instead 
of performing a full inversion. M -1 is therefore not calculated. 
Individual row vectors M^ 1 of M _I can be calculated by solv- 
ing for MMr 1 = where <5, is the Kroneckerdelta. This has 
been carried out, figure [3] shows the results for a part of the 
pixel barrel detector. The sagittae in the local v direction can be 
determined to a precision of a few microns. This procedure de- 
livers the statistical error only. No estimate on systematic errors 
has been carried out yet. 
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Figure 1: Residuals perpendicular to surface along modules: (a) Shown are the observed residuals in the two innermost layers of the strip barrel (tracker inner 
barrel, TIB), expressed as dw = du/ tana. Green circles: alignment assuming a flat surface. Blue squares: assuming a curved surface (2nd order polynomial in u 
and v plus mixed term). The measured quantities were the residual du of the measured hit and the position predicted by the track fit, and the track angle a, measured 
to the normal in direction of u. Only hits fulfilling | tanal > 1/2 have been used. Results from all modules were used, working on 200000 cosmic ray tracks, (b) 
Same shown along v for composite modules in two innermost layers of the tracker barrel with coarser modules (tracker outer barrel, TOB). Green circles: alignment 
assuming one single flat sensor. Blue squares: assuming two flat sensors, splitted at 2v/L v = 0, bows neither determined nor corrected for. Only hits fulfilling 
|tano-| > 1/2 have been used. Results from all modules were used, working on 200000 cosmic ray tracks. The splitted surface assumption results in a flatter 
distribution than the single surface. Observe that in this study there was no correction for the bow, hence the right side shows a bow. 
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Figure 3: Error estimate of some parameters (statistical error only): Shown 24? 
are the sagittae in v of pixel modules in the innermost layer, determined using 
a set of collision tracks (1.3 million minimum bias events at = 7 TeV 248 
corresponding to 2.8 million good tracks) and cosmic ray tracks (2.5 million249 
events, corresponding to 1.8 million tracks selected for alignment). The lower 250 
part shows the error bars centered at zero. 



253 
254 



222 4. Considerations for future detectors 

2! 

223 From the experience of the alignment of a tracking detector 

224 at CMS, the following considerations may be helpful in order 

225 to enhance the alignability of future detectors of similar design. 21 

226 |^] We can only speak for the configuration we know, so these 2 

227 thoughts need proper adjustments for other cases and are far 21 

228 from being universal and exhaustive. 



Resolution. This seems to be trivial. Two aspects are worth 261 
mentioning: The fact that pixels measure two coordinates al- 262 
lows for alignment in all six basic degrees of freedom. Al- 
though strips can be aligned in all three rotations, doing so witli263 
pixels is far easier. A second important thing is the very high 26 4 
precision of the pitch along one sensor ans its constantness on 26 5 
the full area. 26 6 



2 This section reflects the personal suggestions for design considerations for 268 
future detectors of one author (fin) as given in part at the conference. 269 



Module size. Larger size leads to more hits on one module for a 
given spatial hit density. This immediately improves the align- 
ment precision by VN. It also improves the determination of 
angular alignment, as a larger module size translates to a longer 
lever arm. 

Rigid mounting vs. precision mounting. Experiences from 
CMS show, that certain modules can be aligned even though 
they are displaced by few mm (sic!) from the design position. 
If resolution and size are already well chosen, precision mount- 
ing does not necessarily help in improving alignment. There 
might be some configurations for trigger layers, where preci- 
sion mounting may help for other reasons. 

On the other hand, rigid mounting is very important. We 
understand this as that the modules stay in their position over 
time. Track-based alignment needs data gathered over long 
time-periods. As it is averaging in nature, it assumes stability 
during the time required for recording the data it uses. Move- 
ments due to vibrations or imposed by changing conditions like 
temperature or magnetic field must be slower than the typical 
data-taking time. 

Geometric shape. Barrel-and-endcap configurations have a 
great advantage for alignment: They deviate sufficiently from 
an "ideal" sphere-shaped layered detector. The modules are 
also flat (ore just slightly bowed), which naturally leads to a 
spread in incident angles on top of what the event topologies 
may deliver. This helps in creating constraints on several align- 
ment modes. 

Tracks from non-standard origin. Such tracks add more con- 
straints on possible movements of modules which are weakly 
sensitive or even insensitive of changing the x 2 - Cosmic ray 
muons are an example in the case of CMS. They come at large 
do for free whereas in collision data tracks with large do are rare 
events from secondary vertices. They also may connect parts 
of the detector with straight tracks which would normally not 
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270 be connected without imposing special constraints. An exam-321 

271 pie for this are the upper and lower hemispheres of a detector,322 

272 which are connected with a straight track in the case of cosmic 

273 muons. Collision tracks connect these parts as well, but the use 

274 of a common-vertex constraint is necessary. 



thoughts on how future tracking detectors might benefit from 
the experience gathered during our work. 



323 6. Acknowledgment 



275 Optical survey. Survey has a huge drawback. Usually, survey ^ 

276 is performed under certain artificial consitions before the final 

r 326 

277 commissioning. It is a, hopefully well etsablished, assumption,^ 

278 that the survey data stays reliable over time. 

J J 328 

279 Survey is still helpful in several ways. It delivers an indepen- 32g 

280 dent knowledge on the geometry at the beginning of the detec- 

281 tors' operation. As the already mentioned tradeoff is present, 

282 the investment in survey should be limited. Think of an easy 330 

283 way of determining the positions of modules. In our case, po- 

284 sition marks from the layer masks used by the manufacturing. 

285 process were still visible after mounting. A standard single 

286 reflex digital camera with a decent macro lens was used to de 

287 termine relative positions of modules w.r.t their neighbours at a 

288 precision of a few micrometers. Such information can be used. 

289 as an independent measurement for validation od the alignment 

290 or it may be treated as independent measurements included in 

291 the alignment algorithm. 



292 Overlap. Regions with overlap are useful for alignment and343 

293 monitoring of it: Particle tracks have short propagation dis- 344 

294 tances and therefore their trajectories are less prone to effects^ 

295 imposed by multiple Coulomb scattering. The short distance 347 

296 between two sensors along a particle trajectory in reagions of^s 

297 overlap connects them together very tightly. 



298 Unnecessary features. In the case of silicon detectors, imple- 

299 menting hardware-based alignment systems is a difficult task 

300 Either they rely on precision mounting (e.g. some independent 354 

301 sensors mounted on the frame of the silicon sensors) or they 355 

' J 356 

302 mimick tracks by using lasers and holes in the metalization. 

303 Only when their precision is at least comparable to the intrin- 

304 sic track hit resolution, a benefit may be realizable. They also 

305 may suffer from systematic problems, as their tracks have no 

306 geometric spread. 
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307 Alignment studies. The main reason why CMS acheived to 

308 align its inner tracker within that short timeframe was the use of 

309 well-known algorithms, the work of experienced people and the 

310 extensive use of a versatile alignment simulation framework. It 

311 is paramount to have the ability to simulate the detector as close 

312 to reality as possible before the final construction. There will 

313 still be surprises, like the bowed sensors. 



314 5. Conclusions 

315 The inclusion of a more complex surface description of the 

316 slicon sensors of the CMS inner tracker has been shown. This 

317 was able to accommodate for discrepancies found in studies on 

318 the alignment quality and will improve the track reconstruction 

319 in CMS. The sensor bows can be determined with a statisti- 

320 cal precision of a few micrometers. We also presented some 
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